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1 1 3  ABSTRACT 

Developments  resulting  from  work  on  Contract  No.  F30602-71-C-0356  are 
reported.  In  Section  1,  some  remarks  are  presented  on  Furutsu's  recent 
theoretical  analysis  of  irradiance  scintillation.  In  Section  2,  we  define 
and  compute  useful  criteria  for  the  resolution  of  point  features  on  an 
illuminated  target  in  turbulent  air.  The  results  are  related  to  the  work 
on  focal-spot  areas.  The  coefficient  of  the  log-amplitude  variance  in  the 
saturation  regime  is  computed,  and  we  find  <6y^>  =  0.41(Km7'  6L3cn2)-l/6  in 
Section  3.  Finally,  the  work  on  power  spectra  of  angle-of-arrival  fluctuations 
is  expanded  in  Section  4  to  include  a  simple  interferometer.  Two  physical 
cases,  simply  related  to  the  ray  formulas,  appear  to  apply  in  most  practical 
situations . 


DD 


FORM 

*  NOV  03 


1473 


X 


UNCLASSIFIED _ 

Secjnty  Clasii  ‘icatior. 


UNCLASSIF1EP _ 

Security  Classific*ti~m 


1  4 

KEY  WOHOS 


Wave  propagation 
Turbulent  atmosphere,  air 
Laser  beam,  focused 
Optics,  atmospheric 
Irradiance  fluctuations 


The  Quarterly  Report  was  prepared  by  RCA  Laboratories,  Princeton,  New 
Jersey,  under  Contract  No,  F30602-72-C-0486.  It  describes  work  performed  from 
June  19,  1572  to  September  18,  1972  m  the  Communications  Research  Laboratory, 
Dr.  K.  H.  Powers,  Director.  The  principal  investigator  and  project  scientist 
in  Dr.  D.  A.  de  Wolf. 

The  report  was  submitted  by  the  author  on  Id  October  1972.  Submission  of 
the  report  does  not  constitute  Air  Foce  approval  of  the  report's  findings  or 
conclusions.  It  is  submitted  only  for  the  exchange  and  stimulation  of  ideas. 

The  Air  Force  Program  Monitor  is  Lt.  Darryl  P.  Greenwood. 


i 

i 

! 


SUMMARY  AND  ABSTRACT 


Developments  resulting  from  work  on  Contract  No.  F30602-71-C-0356  are 
reported.  In  Section  1,  some  remarks  are  presented  on  Furutsu's  recent 
theoretical  analysis  of  irradiance  scintillation.  In  Section  2,  we  define 
and  compare  useful  criteria  for  the  resolution  of  point  features  on  an  illum¬ 
inated  target  in  turbulent  air.  The  results  are  related  to  the  work  on  focal- 

spot  oreas .  The  coefficient  of  the  log-amplitude  variance  in  the  saturation 

2  7 /6  3  ?  /6 

regime  is  computed,  and  we  find  <  >  =  0.41  L  C^)  in  Section  3. 

Finally,  the  work  on  power  spectra  of  angle-of-arrival  fluctuations  is  ex¬ 
panded  in  Section  4  to  include  a  simple  interferometer.  Two  physical  cases, 
simply  related  to  the  ray  formulas,  appear  to  apply  in  most  practical 
situations. 
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1,  SOME  REMARKS  ON  THE  FURUTSU  THEORY 


In  a  recent  Special  Report  on  the  XVIIth  General  Assembly  of  URSI*  held 
at  Warsaw,  Poland  late  in  August  1972,  we  mentioned  a  discussion  on  Furutsu's 
recent  workfl].  Furutsu,  in  work  that  is  exceptionally  difficult  to  read, 
has  derived  the  following  results  for  beam  waves  propagating  through  turbulent 
air: 

(i)  The  probability  density  of  irradiance  I  approaches  that  of  a 
Rice-Nakagami  distribution  at  large  propagation  distances  L. 

(ii)  Through  another  type  of  approximation  Furutsu  also  derives 

that  the  logarithm  of  I  follows  the  Rice-Nakagami  distribution: 
interestingly  enough  this  distribution  reduces  to  the  log¬ 
normal  one  for  points  far  away  from  the  central  beam  axis. 

In  the  discussion  at  URSI  we  commented  on  these  results,  and  therefore  on 
parts  of  Furutsu's  theory.  Because  his  work  is  hard  to  read,  and  beciuse 
Furutsu's  results  are  being  quoted  indiscriminately,  we  shall  offer  these 
detailed  comments  on  the  theory  as  w^  see  it. 

There  appear  to  be  some  unphysical  features  of  Furutsu's  solution.  The 
above-quoted  result  (ii)  contradicts  result  (i)  for  points  of  observation  on 
the  beam  axis;  they  cannot  both  hold  simultaneously.  Then,  as  Furutsu  makes 
explicitly  clear  in  the  paragraphs  following  hi;.  Rq.  (118),  the  fluctuations 
tend  to  zero  as  the  beam  wave  approaches  a  plane  wave.  It  is  indeed  possible 
that  -  in  some  way  -  the  experimental  realizations  of  plane  waves  are  not  that 
close  to  Furutsu's  plane-wave  limit  of  his  beam  waves,  but  at  any  rate,  it 
seems  decidedly  unphysical  that  the  beam-wave  I  does  fluctuate  on  the  central 
axis,  whereas  the  plane-wave  I  does  not.  We  have  uncovered  a  possible  source 
of  erre’"  in  one  of  Furutsu's  approximations  which  may  make  (ii)  unphysical. 

In  order  to  u'acuss  it,  we  have  to  give  a  very  brief  summary  of  Furutsu's 
work.  It  consists  of  three  distinct  parts. 

l.l  A  Formal  Expression  for  P(I) 

Through  a  method  involving  variational  derivatives  and  cumulants  of  the 
random  field,  Furutsu  derives  a  general  expression  of  the  irradiance  probability 
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density  P(I).  We  quote  a  form  given  by  Bremmer  in  a  paper  presented  at  URSI 
in  Warsaw[2]  for  -plane  waves: 
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a.^  =  Iq  [1  -  exp(-2aL)] 
1/2 

°01  =  ai0  =  Xo  exp(~aL) 


Here.  J  is  a  Bessel  function;  J  (ix)  =  I  (x)  in  real  notation,  but  we  use  the 
o  o  o 

symbol  Iq  for  the  free-space  irradiance  and  therefore  have  utilized  the  Bessel 

function  J  with  imaginary  argument.  The  well-known  formal  attenuation  coeffi- 

0  25/32 

cient  of  the  coherent  (average)  field  is  a  ^  k  L  Lfl  Cn«  The  cumulants 

are  tv.  .  complicated  to  write  down  but  they  are  given  in  both  [1]  and  [2].  The 

summation  over  j  and  k  excludes  j  =  k  =  1.  It  is  this  result  that  yields  (i) 

when  ignoring  all  higher-order  cumulants,  in  which  case  operator  Q(a.,)  w  1. 

3 

Furutsu  and  Bremmer  stress  that  the  off-diagonal  contributions  of  a.,  (for  j  ^  k) 

J 

are  negligible;  this  simplifies  operator  Q(o.,  )  somewhat  but  not  enough  to  add 

J  k 

anything  new. 


1.2  The  Moment  Equations 

Furutsu  therefore  abandons  this  approach  and  derives  the  moment  equations 

by  his  cumulants  +  variational-derivatives  method.  We  state  his  equations  in 

the  following  form:  Let  r  =  (p  ,L)  =  (x  ,y  ,L)  be  the  coordinate  of  a  point  of 
°  m  m  m  'm  „  r 

observation  in  the  plane  z  =  L.  The  N-th  moment < I  >  is  derived  from  a  quantity 

M(pj,p^,L)  =  <?:(r^) - E(r^)E*(rp - E*(r^)>that  is  a  function  of  L  and  of  the 

2N  transverse  coordinates  p^ - ^N'^l - (sy™*50^2^  by  Pj  >  P^)  by  setting 


2 


0,  the 


p'  =  p  ,  and  then  allowing  all  N  coordinates  p  to  coincide  with  p  * 
mm  m  o 

coordinate  of  the  central  axis.  The  moment  equation  is 
fa/3z  -  (i/2k)D  -  (k/2)W(p 


M(p^,p^,z)  «  0, 


(2) 


completely  analogous  to  the  parabolic  wave  equation  in  a  random  dielectric 
where  we  have  the  Laplacian  A  instead  of  operator  D,  and  fie  instead  of  iW. 
The  definitions  of  D  and  W  are 


N 


D  =  iL  (A  -  A’) 


m=l 


m  m 


N  N 

w-  £  L 

j=i  k=i 


v(Vpj)  +  v<pk-pj) "  V(pjA)  ‘  v(p 


(3) 


V(p) 


=  k J  dz<fie(3,z)fie(p,z)> 


Note  that  W  is  wholly  determined  by  the  plane-wave  mutual- coherence  function 

2 

(mcf)  V(p),  which  -  for  the  Kolmogorov  turbulence  spectrum  of  fie  -  goes  as  p 


for  p  <  (microscale) ,  as  p 


5/3 


for  lQ  <  p  <  Lq  (macroscale) ,  and  saturates 


as  p  +  ».  No  objections  to  (2)  and  (3)  and  no  differences  with  other  earlier 
derivations  are  noted. 


1.3  An  Approximate  Solution  of  the  Moment  Equations 

Furutsu  finds  an  approximate  solution  to  (2),  namely  M^(p  ,p^,z)  for 

~y  ~y  -y  -y  -y  ^ 

p^,p^  (all  j  ,k)  by  approximating  W(p_j  ,p^)  by  W^p^p^),  the  function  obtained 

in  (3)  by  using  the  p^  form  in  the  V(p)  for  all  p.  Furutsu's  rationale  for 

so  doing  is  based  on  the  premise  that  because  he  wants  a  solution  to  (2)  for 

the  limiting  case  that  p '■ -p .  -*  0,  for  all  j,k,  he  might  just  as  well  use  a 

*  3 

limiting  form  for  W  obtained  in  the  limit  P^-Pj  0.  That  limiting  form  is 
the  abov< 
in  (ii). 


the  above-described  W  obtained  by  utilizing  Furutsu's  Eq.  (68a),  and  it  results 

r 
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We  believe  the  source  of  possible  error,  or  difference  of  opinion  as  to 

the  validity  of  the  result,  lies  right  here.  In  our  view  the  approximation 

W  ~  Wp  yields  an  Mp  t  M.  To  make  this  plausible,  we  recast  (2)  into  an  integral 

equation  subject  to  the  boundary  condition  that  M  (p.,p’,L )  +  lasW  +  0  (the 

j  K 

f ree-space  moments  are  all  unity) . 


M(P.,%,L) 


=  1  -  (k/2) 


/"[■ 


-(i/2k) (L-z)D 


w(p\.  ,PpM(o^ 


-  ,p£ ,  2 J 


-?j  "  \ 


A  formal  series  solution  to  (4)  is  obtained  through  iteration;  the  re¬ 
sulting  series  is  analogous  to  that  of  the  Born  terms  of  the  integral  wave 
equation. 

Furutsu's  approximation  W  =  W  is  correct  when  D  =  0.  This  can  be  seen 

F  2 

in  general,  but  also  specifically  from  a  simple  example:  Let  W  =  sin  z.  Then 
2 

Wp  =  z  ,  obviously  the  asymptotic  form  of  W  as  z  0.  If  we  solve  (2)  for  this 
simple  case  we  obtain 


M  =  exp  |~(2z  -  sin2z)/4j  ,  Mp  =  exp(-z^/3) 


Clearly,  M  •*  Mp  as  z  -*■  0,  as  Furutsu  claims  it  should,  and  he  was  pre¬ 
sumably  guided  by  such  an  example.  However,  when  D  ^  0,  as  is  the  case  in  (3), 
matters  are  entirely  different.  To  start  out  with,  the  operator  in  (4)  should 

work  first  on  WM;  only  thereaicer  should  one  set  p.  =  "o’ .  The  exponential 

j  k 

operator  containing  D  brings  into  play  an  infinite  number  of  higher-order  deriv¬ 
atives  of  V  at  the  origin  of  its  coordinate.  This  means  that  the  shape  of  V 
(thus  also  of  W)  well  away  from  the  coordinate  origin  plays  a  role  in  determin¬ 
ing  M!  That  suggests  strongly  that  Mp  $  M  because  Wp(p)  #  W(p)  for  p  >  £q[3]. 

We  have,  in  fact,  cor'  meed  ourselves  of  this  fact  by  choosing  a  few  simple 
examples  of  operators  D  in  the  first  Born  term  of  (4). 

Our  conclusion  is  that  Furutsu's  result  (ii)  presumably  is  not  valid  for 
EM  waves  in  turbulent  air.  It  would  hold  onl>  for  a  hypothetical  random  medium 
governed  by  a  V(p)  «  p^  for  all  p.  However,  even  for  such  a  medium  it  .s 
difficult  to  see  how  one  can  reconcile  result  (ii)  with  (i) , 


2.  REf-OLLTIOX  OF  FEATURES  OF  AX  ILLUMINATED  OBJECT 


The  resolution  of  optically  visible  features  of  an  illuminated  object 
in  tht  turbulent  atmosphere  has  been  studied  by  Fried[4,5J.  Fried  has  com¬ 
puted  the  MTF  (Mutual  Transfer  Function)  of  an  optical  system  in  the  turbulent 
atmosphere,  taking  into  account  a  differessce  between  long-term  and  short— tern 
exposure.  His  formalism  describes  the  distortion  of  a  point  image  formed  by 
a  lens  in  turbulent  air.  We  are  interested  in  a  somewhat  different  concept 
and  will  therefore  adopt  another  approach. 

Consider  Fig.  1  where  we  have  sketched  a  ground-level  "eye"  observing  an 


Figure  i.  A  ground-level  observer,  and  two  point  features 
separated  by  distance  pD  from  each  other  on  an 
object  at  range  L  and  height  Z. 


for  small  oq/L.  We  want  to  know  whether  or  not  we  can  distinguish  two  point 

sources  nn  this  object  separated  by  angle  &  in  turbulent  air,  given  that  we 

S 

can  distinguish  between  them  in  tree  space.  To  do  so,  we  shall  compute  the 
difference  63  in  bending  of  tne  two  rays  emanating  from  each  point  source  which 
would  reach  the  eye  in  the  absence  of  turbulence.  This  corresponds  to  Fried's 
short-term  exposure.  Analogous  to  long-term  exposure  would  be  the  distortion 


5 


which  includes  refraction  of  the  central  ray  corresponding  to  inage  blurring 
by  random  variations  of  "he  wavefront  tilt.  This  latter  case  does  not  appear 
very  interesting  to  us  because  it  is  always  possible  to  track  the  total  object 
and  thus  remove  the  '^wavefront  tilt"  as  Fried  terms  it.  The  other  case, 
corresponding  presumably  to  Fried’s  short-term  exposure ,  appears  more  interest¬ 
ing  because  it  tolls  us  hc^  features  be cone  indistinguishable  even  when  the 
total  object  is  being  tracked.  Thus,  we  are  interested  in  the  angle  56  which 
is  appropriately  given  for  all  conceivable  distances  in  the  atmosphere  by  the 
geometrical-optics  ray  formula, 

55  =  —  J  ds  '  Fjeefoj  ( s)  ,s]  -  ^Scio^s)  ,sl^ 

o.(s)  =  (L-s)s  /2L  (6) 

i  o 

p  (s)  =  -  (L-s)p  /2L 
z  o 

The  derivatives  7^  and  are  with  respect  to  the  transverse  coordinates 
p^(s)  ami  p^is)  °*  t**0  rays  (sketched  in  Fig.  1)  we  are  cos aring. 

Clearly  55  -*■  0  as  cq  ■*  0,  and  other  properties  expected  of  this  differ-eme 
angle  58  will  become  apparent  later.  We  shall  proceed  to  calculate  t'-»  sta¬ 
tistics  of  66.  Analogous  to  so  many  other  similar  quantities  in  turbulent  air. 
this  one  too  is  clearly  Gaussian  with  zero  mean,  and  we  need  only  compute  its 
variance.  To  describe  all  the  details  of  this  calculation  would  be  extremely 
tedious  as  it  differs  hardly  from  many  previous  ones.  The  notation  is  given 
in  preceding  reports[6J.  The  technique  is  always  similar.  It  involves  the 
following  steps: 

(i)  Replace  v^ocfpjS)  by  -  iK5t(k,s)  through  a  partial  two- 
dimensional  Fourier  transform  in  a  plane  s  =  const. 

(ii)  Use  the  customary  properties  of  the  covariance <5e(K,s) 

__  ■+  1  i 
oe(K2»S2)>  to  igncre  effects  of  order  L  /L  and  to  introduce 

the  dielectric  turbulence  spectrum  J(K). 


o  7  7  ?  —  1 1  f  o  o 

6  ;  (K)  =  32r  X  0.033  C  (K  +  L  )  exp  (-K  /O 

no  m 


(7) 
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We  f lad 


L  * 

<®32>=  (4x)_1J^  dSf 2(s) d£  K3  4(X)  Jl  -  J^O^Cs) |J 


(8) 


5 

l 

s 

i 


where  p  (s)  is  given  in  (6).  As  we  shall  consider  fairly  small  values  of 
^  -> 


p. (s)  and  because  there  is  a  factor 


fa  the  above  integrand,  it  can  be 


A  2  — H/3 

argued  that  the  factor  (X  +  L  )  in  (7)  nay  be  replaced  by  K  .  Its- 

°  2  2 
serting  (6)  and  (7)  into  (8),  introducing  the  auxiliary  variable  x  =  K  /<  , 

a* 

and  then  performing  the  integration  ever  the  new  variable  x,  we  obtain: 

<S6^>=  1.3  r(i/6)«^/3  J*  ds  C2(s)jl  -  m[i/6,1,-<2c2(L-s>2/4L?]|  (9) 


The  coefficient  1-3  is  an  approximation  of  4  *  0 . 03 3t  accurate  to  about  0.22, 
F  is  a  gamma  function,  and  M(=,b,z)  is  a  Kramer  function  defined  in  (13-1.2) 
of{7]. 

2.1  Horizontal  Propagation 

2 

tor  horizontal  propagation,  Cq(s)  is  a  constant.  Using  the  above-cited 
power-series  definition  of  the  Kuaraer  function,  we  expand  into  a  power  series 
and  integrate  tern  by  tern  to  obta in 


<552>  =  1.3  C2e1/3L  x  £ 
11  13  13=1 


(-D^rq/bfrB) 

(n!)"(2irH) 


■(*) 


2a 


(10) 


This  fora  is  useful  for  snail  values  of  <  P  .  ihus,  when  p  «  f.  ,  i.e.  when 

no  o  o 

two  distinct  features  are  separated  by  much  less  than  the  Bicroscale,  we  find 
a  differential  refraction  with  variance. 


<c~Z>  =  1.3  cV/3L  x  1  (<  c  )2  for  p  «  i 

n  m  12  b  o  o  o 


(11a) 


Note  that  the  coBbination  i  3  v.2*c3/^L  is  identical  to  the  ratio  r2  /L2  where 
2  n  b  LB 

-r  is  the  area  or  the  average  focal  spot  in  turbulent  air  minus  the  vacuum 
LB 
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diffraction-limi <ed  area  as  defined  m  RAbC-TS-72-il9[6}-  On  the  ocher  hand, 

vh&n  p  »  1  the  Sumner  fuacticn  in  (9)  is  negligible  and  we  obtain 
o  o 

<552>  =  1.3  C2jex?3L  *  F(l/6)  =  r(l/6)rfa/L2  for  a  »  1  (lib) 

a  he  o  o 


Equations  (11;  provide  a  basis  for  deciding  whether  or  not  revolt  cion  is 

2^ 

possible-  vs  suggest  -  somewhat  arbitrarily  -  that  <c-J  be  compared  with 
*  2  2  2 

S  '-q/L.  The  reason  for  this  criterion  is  that  a  ray  coating  from  feature  1 

go 

would  appear  to  be  coming  from  feature  2  when  these  angles  are  equal,  aad  hence 

feature  2  would  not  be  indistinguishable  from  feature  1.  However,  features  1 

7  2 

and  2  are  not  confused  with  each  other  when <©S  ^  «  &  ,  i.e.,  resolution  is 

2  g 

good.  A  resolution  index  of  use  would  be  <(55/5  )  >  .  Resolution  is  good 

S 

(bad)  wnen  this  index  is  much  less  (more)  than  unity.  From  (11)  we  obtain 


<(5^/5  )2>  = r  )2  for  a  «  1 
g  12  n  Lr  o  o 


=  r(’/6)  (r  /p  )  for  p  »  1 
LB  O  O  o 


Consider  (12)  first  of  all  for  stjcII  p  (p  «  i  ),  i.e.  for  resolving  well 

o  o  o 

within  the  microscale.  It  is  noteworthy  that  the  criterion  is  then  indepen¬ 
dent  of  pq!  In  fact  it  simply  depends  upon  the  ratio  of  r^g  to  i^,  i.e.  upon 
the  strength  of  cumulative  turbulence  along  the  path,  whether  or  not  resolution 
of  small  features  is  possible.  Let  ns  recast  (12)  into  a  more  numerical  form: 


<(6e/e  )2>  *  (r/o.6i  i  )l 


i  )  >  -  (r7_/U.bl  l  )  r or  c  <<  t 

g  L3  o  o  o 


*  (2'4rLB/Co)' 


for  o  »  i 
o  o 


We  have  reproduced  Fig.  1  of  RADC-TR-72-119  as  Fig.  2  in  this  report,  so  that 

2 

the  user  can  simply  read  off  r.  _  for  given  L  and  C  (l  was  chosen  to  be  6  mm 

LS  n  o 

for  this  graph).  Sub-microscale  details  appear  to  be  unresolvable  for 
2  -14  -2/3 

C  >  30  m  and  distances  beycnd  several  hundred  meters,  a  hardly  surprising 


Figure  2.  A  plot  of  ryj  as  a  function  of  range 
L  and  turbulence  strength  C^. 


fact  and  somewhat  useless  since  the  microscale  is  in  itself  so  small.  Let  us 

rather  consider  the  other  -  more  practical  -  case  that  c  »  h  in  which  case 

o  o 

the  second  of  (13)  would  hold.  In  this  case  resolution  is  i->st  for  all  details 

that  exceed  the  size  2.4r  which  is  easily  read  off  Fig.  2.  For  example,  the 

Ld  2  2/3 

critical  size  for  distinguishing  details  at  3-km  distance,  when  c  =  10 

n 

appears  to  be  of  the  order  of  8  cm. 

A  more  formal  version  of  this  simple  theory  can  be  obtained  by  utilizing 
Yura's  procedure[8]  for  cor  uting  t:.e  mutual -coherence  function  (mcf)  of  two 
deita-functior.  sources  at  locations  1  and  2  of  Fig.  1.  It  appears  sufficient 
to  us,  however,  to  utilize  (13  and  Fig.  2  for  rapid  estimates  of  the  amount 
of  .etail  that  can  be  resolved.  In  the  sense  that  amplitude  effects  do  not 
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appear  to  play  an  important  role  in  the  plane-wave  and  spherical-wave  mcf,  it 
would  appear  that  considerations  based  on  (6)  suffice. 


2-2  Slant-Path  Propagation 

If  the  viewed  object  is  at  altitude  2,  then  the  development  from  (8)  cm 
2  2 

is  different  because  z  (s) ,  hence  C  (s),  is  a  function  of  path  position.  The 
2  a 

altitude  model  of  C  (s)  is  discussed  in  73-ELAi5C-72-'ii5f  6]  in  enough  detail. 

n 

ife  shall  adant  the  s  »  l  .audition  so  that  the  Sessel  function  in  (8)  can 

o  o 

be  ignored.  In  that  c~se  we  find 
2 (XT'  =  ?(Z>  <c«2(0)> 

1  (14) 

?W  i  as  [  (1-5-  sZ/z  )"2/3  (1  -s-  sZ/1  >_2/3  ezo(-2sZ/h), 


where  '4S  ('<)>  is  the  horizontal-proi  agation  result  given  L>  (lib).  The  Xonin- 

Cbukhcv  ien-  th  has  been  divided  by  7  to  yield  1  ;  at  dawn  and  dusk  I  -»■  =  and 

s  s 

at  midday  (sunny  weather)  I  *1.5  i.  The  observation  a>titude  is  chcsen 

s 

arbitrarily  as  z  ==  1  a,  but  there  is  some  uncertainty  ia  this  choice.  There 
o 

is  also  considerable  arbitrariness  in  the  use  of  this  model  for  altitudes 
above  several  hundred  meters.  As  a  result,  we  shall  also  restrict  ourselves 
to  Z  «  h. 


(A)  Dawn,  dusk  approximation,  £.  -»•» 

In  this  case,  only  the  first  integrand  factor  of  F(Z)  deviates  from  unity, 
,.nd  we  obtain 


■(Z)  --  3  (Z/z  )~z/i ,  ?  «  Z  «  h 

o  o 


(15a) 


(B)  Kiddav,  sunshine,  >.  =■  1.3  m 

s 

This  is  a  slightly  more  complicated  case,  even  for  z  :7  <<  h. 

c 

coordinate  t»  ctnsf  ormati ons ,  C14)  can  be  recast  into 

ZK.ls-z0) 

T.2  2  3  .  1 1/3  f.  -2/3  ...  -2/3 

F(Z)  LVc'22  *' VZo)J  J  dy  y  U+>'-* 

0 


Through 
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which  expression  contains  cn  integral  representation  of  ia&wplete  beta  function. 
Ecwever,  toe  upper  bound  of  the  integral  is  extremely  large;  hence,  we  way  set 
it  to  <*.  We  then  have  rbe.  representation  tf  d«  i'saaner  fsocticn  G(l/3,  2/3,  0) 
a  (sec  formula  »'I3.2.5)  in  ref.  f 7 j> ,  afeich  restates  to  a  quotient  of  products  of 

f  gamma  functions.  We  flag 


r(a  ^rq/j) _ 

'Z'  F(5/6)  .1/1. 


,2/3 


.-1 


z  "  {£  -2  ) 
o  so 


ITT *  c,‘.)  ■  *„ 


«  Z  «  is 


^  5(Z/zc) 


-i 


(15b) 


The  numerical  estimate  in  (15b)  is  bases!  on  l  <1.5a,  2  *1  n.  These  formulas 

s  o 

(15)  are  intended  to  be  guidelines  rather  than  a  rigorous  basis.  Toe  uncertaia- 
2 

ties  la  the  C^(s)  model  are  too  large  ,or  more  accurate  results.  We  note  a 
stronger  reduction  with  altitude  increase  of  <36i(Z)>  Caspar ad  with  an  equiva¬ 
lent  horizontal  -path  propagation  situation  at  midd  ty  than  at  dawn  or  dusk;  on 
2 

the  other  hand,  0^(0)  is  usjaiiy  mich  higher  «;t  midday  so  that  conditions  for 
feature  resolution  are  usually  better  at  dawn  and  dusk,  even  for  objects  fairly 
high  in  elevation  angle. 


11 


3-  COEFFICIENT  0?  THE  L0C-*SGpLiT53E  TASIASCE 


A  new  theory  of  irraaicnce  fluctuaticcs  of  a  plane  wave  la  turbulent  air 

was  presented  la  a  previous  quarterly  report,  £A3C-TlE-7?-204f 6]-  Its  major  new 

resale  was  to  predict  a  log-normal  distriSratloa  of  irradiarce  1  la  the  satura- 

tloa  regime,  l.e.  where  k?  is  a at  small.  Xiore  specifically ,  it  de— 

2  ?  j/-j  ?  2 

fiaes  true  saturation  regime  by  the  size  of  the  parameter  j‘  =  <•'  kl.  C^,  which 
does  ant  differ  a  great  deal  freer  ^C2  at  cotical  frequencies,  aad  it 

ES 

predicts  for  the  log- amplitude  variance 


<ii2>'  « 

m  a 


(17) 


when  rT  »  1.  The  right-hand  side  of  (17)  may  Icok  unfamiliar  but  it  too  is 

mot  very  different  front  */o  at  optical  frequencies.  One  should 

r  ’  2 

note,  however,  that  the  predicted  comcinatxons  of  variables  is  aad  la  (17) 

are  not  restricted  to  Kolmogorov  turbulence  statistics.  They  are  more  general , 

and  st  prefer  to  use  them  also  because  the  distinct  ion  between  them  and 
.7/6,11/6.2  .  . 

at  L  C  cannot  oe  made  experimental  Iv  to  cate. 

n  *  ' 

A  coefficient  of  proportionality  is  missing  in  (17).  The  purpos »  of  this 
section  is  to  estimate  it.  Toe  starting  point  is  Eq.  (A4)  in  5ADC-TH-72-204- 
This  equation  yields  a  real  part  x  giver,  by 


X  =  -(k/8-r2)  J"  dz  j"  dZ  (K,  z)sin{K“(I.-r)/2k]exp(-iK*c  (z)  ] ,  (18a) 


and  we  use  the  conservation-or-energy  relationship  for  the  variance  of  y , 

2  2  _ 
namely  <  iy  >=  _  <>.  >,  to  compute  the  variance  <c.v  >  iron  (18a).  the 

correlation<d&(K,z)exp[-iK*o(z)  ]  >can  be  computed  to  oe 

(eV4)K2J (K)(L-z)exp[-K2<c2(z)>  /4]d2K.  Substituting  this  into  (13a), 

2-11/3  22 

utilizing  the  spectrum  *(K)  =  32yC  k  exp(-K  /<  and  making  the  variable 

2  n  h  ^ 

changes  of  K  to  r.  =  K  (L-z)/k  and  z  to  y  =  (L-z)/L  we  obtain 


<CX2>  =  |'fk7/6Ln/V  J  &7  75/6  f  dx  xmex?l~xs*Hy)}, 


(18b) 


where  y  *  O.Q33t  =*  0-325  a  rad  f(y)  *  yr(7/6)y(3-2y) .  The  last  integral  ia 
(18b)  yields  a  gams  fooctica  F(7 /6),  and  upon  sts«*  assent  rearranging  of  the 
factors  we  obtain 


<^>-4J'W) 


-  r  ■*  f-t iej\ i”7-/ 6  I  f  .  -1/3,,  »  «-7/6 

£  *  Ifi (7/6)3  *2  7  fly  7  (3~2y) 


The  coefficient  2  is  easily  computed  numerically,  and  we  finally  obtain 


.2^  n  *,/  7/3_3_2\-l/6 

>  •  ‘-“I',  *-  CJ  • 


This  is  the  final  result-  In  order  to  facilitate  comparison  with  existing 
data  plotted  as  measured  <©y^>  vs.  observed  values  of  -  0-307k^^L^^C^ 

c  U 

ve  rewrite  (20)  as 


2  0.34(A/k)  7/36 (oe2)' 


This  formula  indicates  a  t/go*  dependence  a  frequency  and  pathlength  of  the 
coefficient  of  (c?)  If  we  set  1  *  6  mm  (i.e.  <n  ~  10^a  ^),  X  =  0.6  ym 

(i.e.  k  =  10^a  *),  and  choose  T  =  1  km,  we  find  the  coefficient  of  (o^) 
to  be  014.  If  we  take  L  =  500  m  or  set  >.  =  0.3  urn,  we  obtain  0.16. 

2 

Mevers,  et  al.[9]  present  two  curves  in  the  saturation  regime  of  <ox  >  vs. 

O 

"7.  In  our  notation  they  find 


■>  •>  -o  17 

<fx“>  =  0.55  (s“)  U*i/ 


at  X  =  0.514  um 


_  ,2.  -0.17  at  X  =  1.15  urn 

=  0.24  (c{) 
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The  agreement  of  the  theory  with  these  dati  is  not  unreasonable,  but  the 
predicted  wavelength  dependence  of  the  coeffl ;ient  is  strongly  contradicted. 
However,  the  1-1 5— ten  curve  is  based  on  only  one  point  in  the  saturation  regise 
so  that  perhaps  one  cannot  use  these  data  to  predict  the  wavelength  dependence 
of  the  coefficient.  More  recent  data  taken  by  Kerr[10]  -•  see  his  Fig.  7  -  in¬ 
dicate  a  weak  wavelength  dependence  hidden  by  data  spread,  and  a  relationship 

<®X^>  ~  0.3  (©^)  (error  perhaps  30%),  (23) 

where  =  0.124  £s  the  parameter  for  spherical  waves.  We  have 

not  yet  computed  our  coefficient  for  spherical  waves  so  that  we  have  net 
checked  out  this  coefficient  which  contains  an  appreciable  error  due  to  data 
spread.  However,  it  appears  to  be  the  sane  order  of  magnitude  as  for  plane- 
wave  data  -  see  Fig.  3  of  Gracheva,  et  al.[LL].  For  these  latter  data  we  infer 
from  the  plane-wave  plot, 

<6x^>  ~0.25  (o')  (error  parhaps  20%)  (24) 

so  that,  again,  there  appears  to  be  reasonable  agreement  with  our  results  (20) 
or  (21). 
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4.  ANGLE-OP-ARRIVAL  POWER  SPECTRUM  FOR  AN  INTERFEROMETER 


In  a  previous  report  (RADC-TR-72-119) [6] ,  the  power  spectrum  of  the  focal- 
spot  scintillation  was  shown  to  be  that  of  the  angle  of  arrival  of  a  ray,  and 
it  was  computed.  In  the  present  report  we  shall  do  the  same  for  an  interfero¬ 
meter  simplified  to  its  essentials.  Basically,  our  interferometer  consists 
of  two  parallel  rays  separated  by  a  distance  p.  They  are  collected  by  a  lens 
at  distance  L  and  focused  into  a  fringe  pattern.  The  "center"  of  this  fringe 
pattern  deviates  from  the  central  axis  only  through  mutual  ray  bending  caused 
by  atmospheric  irregularities  much  larger  than  p  that  are  intersected  by  both 
rays.  All  other  irregularities  give  rise  to  relative  bending,  i.e.,  to  dis¬ 
tortion  of  the  fringe  pattern.  A  measure  of  the  distortion  is  given  by  the 
difference  60  of  the  two  angles  at  which  both  rays  arrive  at  the  lens.  This 
difference  is  described  by  the  vector. 


St 


L 

(1/2)  /  dzf^T6e(p/2,z)  -  ^x6e(-p/2,z)] 
°  L 

-(i/4Ti^)  J'  dz ^*d^KKog(K,z)  sin[K*p/2] 

0 


(25) 


Each  ^  refers  to  the  gradient  operator  working  only  on  the  first  two 
transverse  coordinates  in  6e,  denoted  by  +p/2.  A  partial  Fourier  represents- 
tion  has  been  introduced  in  the  second  line  of  (25).  The  vector  p  describes 
magnitude  and  direction  of  the.  separation  in  the  plane  of  the  two  rays.  Time 
t  is  not  indicated  explicitly  in  6$  or  in  6e(K,z). 

Let  us  define  <J^(K;r)  as  the  three-dimensional  Fourier  transform  with 
respect  to  r^-^  °f  <Se(rx,t)6e*(r2,t+r)>  .  Let  us  also  define  the  auto¬ 
covariance  Rq(t)  =  <6(f(t) *6$*(t+T)  >  .  Upon  forming  this  from  (25),  making 
the  customary  approximations  for  L^  <<  L,  and  performing  the  usual  other 
manipulations,  we  find 


R0(x)  =  (e^L/8n^)  ^  [l-cos(i£*p)  ] 


(26) 
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Note  that  this  expression  differs  from  Eq.  (28)  of  RADC-TR-72-119[6]  only 
by  an  extra  factor  which  describes  the  interferometer  situation.  We  now  apply 
Favre’s  hypothesis  for  $^(K;t), 

$4(K;t)  «  $(K)  exp[iK-?TT  -  4(KAUx)2/3]>  (27) 

as  set  forth  in  RADC-TR-72-119[6] ,  noting  that 

e2$(K)  =  32ttyC2(K2  +  L“2)"11/6exp(-K2/K2)  (28) 

Y  5  0.  033tt2  «  0.325 


We  insert  (28)  and  (27)  into  (26)  and  take  due  note  of  the  fact  that 
there  are  two  fundamental  vectors  in  the  problem:  (i)  the  transverse  velocity 
vector  U^,  (ii)  the  transverse  separation  vector  p.  Before  writing  down  the 
result,  we  develop  notation  some  more.  Let  X  =  p/Lq  be  a  normalized  separa¬ 
tion  vector  (i.e.  in  terms  of  the  number  of  macroscales).  Let  £//  be  its 
component  in  the  direction  of  iL,  and  be  the  component  perpendicular  to 

U  .  Let  us  also  introduce  the  two  angular  frequencies  u  =  n/L  and 
X  X  T  G 

Aco  =  4AU/L  The  result  is  then 

o 


1 
(29) 


00 

R0(t)  =  4yC^LLq  j  dx  x(l  +  x)  11^exp  | -x  j^(Au)x/2)  2+((c^Lo)  ^jf 

0 

V  f”  T  ,  1/2,  1  T  ,  1/2,  1  .  ,  1/2, 

X  [VV*  }  “  2  Jo(p+  x  >  "  2  J0(P~  X  > 

p±  fv  ±  V2 +  ?i2] 


We  now  define  the  power  spectrum  W0(w)  as  the  integral  from  t  =  0  to  x  =  » 
of  2R0 (x) coswx  as  in  Eq.  (32)  of  RALC-TR-72-119[6] .  This  yields 
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1 


w0(«)  =  8ycjLLo'1/3tW0(a),5/6)  -  W0(u>,ll/6>] 


00  00  2 
WQ(td,q)  =  .-1/  dy  cos  y  j  dx(  1  +  x)"qexpj  -x  +(lcmL0)~^]| 


(V1)  - 1  Jo  (?V1/2)  - 1  ^  &-*m) 


21  1/2 


y+  =  [  (y±?/  u/ojt)2  +  (Ciw/a)T)2]  ] 


This  form  for  W  (u>,q)  can  be  simplified  somewhat  by  a  coordinate  trans- 

O 

formation  y  -*■  -y  for  the  third  term.  After  some  algebraic  work,  we  obtain 

00  00  ^ 

W0(w,q)  =  a)  1  J  dycosyj'  dx(l  +  x)  qexp  j  ~x[(~  *  +^lcinLo^~2  | 

°  °  (31) 

X  |  Jo  Gp™172)  '  -Os(^a)/a.T)Jo  (^[y2+(q  a)/a,T)2]1/2x1/2)| 


Note  that  the  parameters  =  Pj^u/U^,  and  £,//  co/o>T  =  to/UT  are  the 

products  of  to  with  a  fundamental  time  T  ^  =  P^/UT  (and  =  p^/U^,)  that 
describes  the  time  required  for  air  to  flow  across  the  interferometer  beams. 

A  third  result  ensues  for  the  case  that  £i  *  0,  i.e.  for  the  case  that 
U.j,//p.  The  obtained  spectrum  is  identical  to  that  for  a  ray  multiplied  by 


tne  factor  [1  -  cos(<oT//)]  where  T^ 

.  i  _  r i _ .  i  t  _  c _  m^l 


' /U^. .  It  filters  out  all  harmonics  of 


the  fundamental  beat  frequency  T^  11  U^/p^. 

From  here  on,  we  will  set  £//  =  0  to  obtain  less  trivial  solutions  for 

air  flow  in  chat  plane  through  the  central  axis  that  is  also  perpendicular  to 
the  plane  of  both  interferometer  rays.  It  is  quite  easy  to  amend  the  solutions 
thus  obtained  to  the  general  case,  simply  by  prefixing  the  factor  cos(toT^)  in 
front  of  the  result  for  the  second  term  of  (31). 
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4.1  Frozen  Flow  (Taylor's  Hypothesis) 


Let  us  first  solve  (31)  for  the  case  that  AD  =  0,  i.e.  no  randon  coapcnent 
of  velocity.  We  obtain. 


OO  CD 

W0(o),q)  =  a  1  J°  dy  cos  yj*  dx(l  +  x)  q  erpt-xGc^)  *] 


-  Jo('«r[5'2  +  a<“'“  t)'i1,2’1/2)| 


(32) 


X  <J 


Following  Appendix  B  of  Clifford's  work  for  the  y  integration[12] ,  we 
reduce  (32)  to 


co 

w0(“»q)  =  f  dx(l  +  x)  q  (xoj2  -  (O2)172  exp (-x/ k^L^) 


(co/o)T)  ■ 


jl  -  cos(5i  [x  -  (w/(dT)2]172)  j 


(33) 


The  following  steps  are  carried  out  on  (33) :  First  substitute 
2  2  1/2 

x  =  z  +  a)  /ojx.  Then  expand  cos(Lz  )  into  its  power  series,  substitute 
2  2  ^ 

z  =  t(l+co  /o>T) ,  and  apply  formula  (13.2.5)  of  ref.  [7]  to  the  integrals  in  t. 
The  ultimate  result  is 


WQ (m,q)  (l  +  a)2/(o2)1/2“q 


-u2isi 

e  T  U  I  1/2,3/2-q,  --r 

1  a: 


2,2' 

m  +  a>T 


[  CO 

>1  -  T  JL 

'  m=0  m! 


r  2 

(i  + 

4  r  2 

“t 


m 


(34) 

2  .  2.  #„2, 


U £  1/2  +  m,3/2  -  q  +  m,(u>  +  a 


y  [3/2, 3/2,  (to2  +  to2)/f!2) 


quite  clearly  an  extension  of  our  previous  result,  Eq.  t36a)  in  RADC-TR-72-119[6] . 
In  the  general  case  that  £,//  /  0,  the  summation  sign  in  (34)  will  be  preceded  by 
the  factor  cos(£^  oj/'u)^)  .  The  Kuramer  function  U(a,b,z)  is  given  in  formula 
(13.1.33)  of  ref. [7],  The  power  spectrum  has  three  distinct  regimes.  For 


IV'.SIJI 


“  '<  i>>,  it  is  flat.  For  oil.  «  at1  «  SL  there  is  a  decrease  in  b\(<u,q)  wits 

1  O  *  •  ^J1  ** 

(“/“V  q  (note  that  =  'mV^*  Fro®  **“  first  of  Ej?5-  (30)  we  note  that 

scon  after  a s  increases  above  the  c  =  5/6  term  is  dcmlaaat  because  it  y  <lfe 

1  -2/3 

a  power-spectral  decrease  as  (ar/o^)  which  certainly  masks  the  small  correc¬ 
tion  from  the  q  =  11/6  term  yielding  a  (ss/as^)  d  crease.  Thus ,  the  sower 
spectrum  decreases  as  as  ^  in  the  regime  of  greatest  interest; 

°X/Lo  <<  “  <<:  °t/£o*  When  “  exceeds  i-e.,  when  we  look  for  spectral 
components  determined  by  eddies  of  smaller  size  than  the  microscale,  we  ob¬ 
serve  a  rapid  cutoff. 

With  regard  to  the  interferometer  "correction  factor”  { - }  in  (34) , 

it  is  obvious  that  W 0(m)  -♦■  0  as  Ci  4  0.  It  is  less  clear  what  happens  when 
^1  00 *  know  that  the  interferometer  separation  becomes  infinite,  and 

because  eddies  are  blown  across  the  interferometer  in  planes  perpendicular  to 
that  of  the  interferometer  (which  are  parallel  to  the  cctlral  axis)  we  know 
that  independent  realizations  intersect  each  axis.  Consequently,  the  factor 
{  }■*■!>  and  we  obtain  twice  the  single-ray  power  spectrum.  This  can  be 

shown  more  easily  by  firs:;  letting  »  in  (34),  i.e.  by  considering  the 

asymptote  of  o>  «  In  that  case  (34)  reduces  to 


WQ(a),q) 


T(q  -  1/2)  /T 

T(q)  u>T 


\ 

) 


1  - 


23/?-q 

r(q  -  1/2) 


cos(^u/a-T)^iq  1^2 


q/2-1/4 

Kl/2-q 


(35) 


It  is  quite  clear  from  the  asymptotic  form  of  K  as  5,  «  that  { - }  -s-  1. 

Furthermore,  (35)  is  clearly  much  handier  for  numerical  purposes  than  (34) , 
which  nas  more  detail  in  it  than  we  really  require  in  the  interesting  portion 
of  the  power  spectrum  (<o  <  Q^) . 

4.2  Random  Flow  (Favre’s  Hypothesis) 

Consider  the  general  case,  AU  *  0  and  LT  /  0:  the  case  of  random  Gaussian 
velocities  with  a  mean  U^,  normal  to  the  propagation  direction.  We  have  not 
been  able  to  perform  both  integrations  in  (31).  Possibly  one  of  the  integrations 
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i^iW^fl8<l.W)WW>i,tyi(Al..f^WtyiWiWM.H'Sy*WWl?l»@88WRWWC8>IHSH!IS!H55?^WI^fl^(W 


c*2  be  c*rriei  eoc  witherac  simplifying  che  expressiccs.  E&werez,  ft  is  possible 
Co  obcafa  same  c loseS-fosa  resales  ey  ahasdsaieg  of  irfry  Satcef Is.  Sfe  vail 

discoss  s  few  of  s£<*se. 

—2/3 

*4.2.1  Use  <u>  Portion  of  rbe  Spectre® 

In  Cfeis  sabsectica,  we  restrict  oorseives  to  ai^  <<  »i  <<  2^.  or  iu  <<  «  «<  £& 
0&2  is  defined  sianrlariy  as  dm  with  *SF  replacing  E—}.  This  rr b&ir  we  set 

U  f  z)  •  *  x  ,  thee  sec  q  -  5/6  and  ignore  else  o  =  11/6  coacrioatiioa,  and  let 
cm  -*■  *  Is*  (31)-  5fe  change  variable  a  to  z  =  t  ,  **»«  obtain 


S»s(«,5/6)  =  1(0)  -  cosUfto/u^JliZ,} , 

*  «t 

I(5I>  =  2m  1 J  dy  cos  y^  dt  t'2/3Jojt[y2  +  (Cs  w/a^) 2]i/2  u^/t  * 


(36) 


x  erpl-t2(y-5i£/2is)2] 


The  following  algebraic  steps  are  performed  upon  (36) :  (i)  multiply  the 

integrand  by  1  =  exp[  (C/^/2i3^,)2t2J  ezp[-(£j  &is/2s£^.)2t23  and  absorb  the  last 
exponential  into  the  last  factor  of  (36);  (ii)  expand  expU/  ks/is^t2)  into 
its  power  series;  (iii)  perform  the  dt  integration,  using  formula  (11.4.28)  of 
ref. [7];  and  (iv)  perform  the  tena-fer-tena  dy  integration  using  formula  (9.6.25) 
°f  ref. [7J.  The  result  is 


1^i)  ~  u  ml  (4t}2C1/‘«>Tb»)1/3(Ci«/2<£T)ID 


x  H(l/6  +  mJl,~i4/6U)2)Km_1/3aio)/ti;T) 


(37) 


where  the  Kummer  function  M(a,b,z).  rlso  known  as  1F]L(a,b,z),  is  defined  in 

(13.1.2)  of  ref. [7],  and  K  (z)  is  a  modified  Bessel  function  of  fractional 

v  —2/3 

order  v.  The  power  spectrum  in  the  w  '  region  is  thus  in  genera), 
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r 


».  Cm)  *  STfCTJL 


C-i*,5/6) 


\  -2/3 


2  ? 


V^,5/6>  *  rCl/3)  SCI/6>,I,-4/£w~> 


j  2 sasUftaW  f*  l_  ('££l.\ 

r  v*J 


aefrI/3 


SoCe  chat  the  first  farter  of  St._(i»,5/6>)  is  in  agreement  with  (vi),  7* 

■*  -» 

of  2A2BC—TS-72— 2©£f6j  ;  there  is  only  a  difference  in  noma! i ration  by  a  constant 

7  — i  ?  —a  ^3 

factor  SvC  ii_  '  (the  factor  4-yC  id.  ~f  J  was  inadvertently  emitted  in  the 
s  o  no 

ray  formnla) . 

Tavlor’c  hypothesis  is  obtained  frem  (38)  in  the  limit  dm  •*■  0.  It  can 

2  2. 

be  seen  that  (35)  resales  vita  only  one  difference:  the  factor  1  -r  u  /»_  in 

2  2 

(35)  is  replaced  by  ii  /u,.- 

Ihe  other  limit  is  -»  0.  However ,  oar  formulation  leading  to  (38)  is 

really  not  suitable  for/ dealing  vita  this  limit-  As  s^.  -  0,  there  is  no  longer 

a  preferred  direction  along  which  a  £/  and  C i  caa  be  defined.  The  lector 

cos(-;  /  /  )  in  (33)  appears  moreover  to  become  infinitely  oscillatory  and 

hence  acts  effectively  as  a  zero  factor  -  unless  C  /  0  more  rapidly  than 

s-  -»•  0.  Clearlv,  the  formulation  in  this  fashion  is  not  suitable  for  small 
I 

ratios  of  _/Ai  or  of  The  difficulty  lies  in  the  transition  from  (30) 

to  (31):  the  formula  where  the  cos(C /  ts/ti-)  factor  first  appears. 

Reconsider  (30)  and  rewrite  the  argument  of  the  second  and  third  Besse. 
function  terms  of  the  integrand  as 


(39) 


I&ebs,  JLai  this  Bfaciging  case,  we  ffiao  irttie:  cfee  C33~ly>-sra*ricggad  ffart sz  off 

j  —  t>  f  »j  2  2 

(2®5  redhxras  8®  1  -  J  C£«sJ  ■*•  ©Co^j/aOJ  »  arn£,  iff  ©Caiyriu  >  casas  are  ignored, 

a*  c&serwe  gfeat  frsotad  off  (2&5  we  arw  occaia  after  cfee  dr  irregraEffotr: 


«sfe,5/6)  «  !•(©  -  I'CO. 


I"  CO 


u  \  2ui 


£_  l»\  -m  f  £z 


-2/3  -c 


J 


£  J  (£u/£tu 

O' 


s*^2) 


*4€a)! 


Ms  ffaunmaia  cant  ce  fcrtber  developed  imt©  a  power  series  off  C^/lui)-,  bet 
tbere  is  no  aeiaacage  t©  so  for  gemerffcai  wcri-  ttfe  simply  compete 

I*  (©)  to  ©bears 


3L(u,5/6)  =  "(1/3)  (/r/&u)  (i:/ixa)'2/j> 

*  {»  -  ni75)  /  i--  =1/2>; 

0 


(4Sfo) 


as  the  leading  term  la  tae  small  iSy/iat  Unit  wiaea  u  »  £w. 

4.2.2  Tfee  Fiat  Portion  of  the  Spectrum 

Here,  we  restrict  ourselves  to  u-  «  !xay  is>  «  u^.-  ite  set  -*•  =,  replace 
y  by  variable  y 1  via  y  =  y*  (then  »e  drop  the  prime  for  convenience ;  and 

then  let  u  -»■  0-  Ve  obtain 


W=(0,q)  =  1(0)  -  cosCCjte/iOKg:), 

v  I  * 


(41) 


Kd)  -  -T_1  f  d,  y  did  +  *)-■>  ^y2  + 


Unfortunately,  we  are  not  able  to  reduce  this  expression  greatly.  By  using 
the  addition  theorem 
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j  a?2  *  V/2-*>  *  j  a  «4^ 

O'  •  tS?  ~  ®  aK  •  £K- 


cast  recast  (41)  faffig  th*  f<s*x 

/"  ~uilJ2J!L^  *  7  7 

I«  )  «  f1  «  ;I  «4/2&OF  CC:  .0 

*■  J  O'  !!  (S'  *■ 


*2E  i-D^ii^/ZLa^WJi  L,<$) j 

»*i 


(43) 


FJ 


?_,«>  5  J"  &E  **£/2Cl  * 


THoe  fTasaccigsas  I  la  f43*  are  modified  iessal  facactinras.  The  fsactiiOSK  F  (£  ,c) 
a  st  ~  " 

«s*r ;  £0  be  ©amperee  saaericallT..  Sowerer,  1(0)  is  easily  cswpsted  because 

F  (3,®)  =  0  for  at  #  0  are  F  (0,c>  =  T(c  -  l/2)/r(e).  Ue  dtcs  ftaf  for  rise 
a  o'" 

flat  part  of  the  spectra®  when  *  0 


¥_(Q)  =  s-fchi  1/3  x  cr(i/3)/:(5/6)  -  r{4/3)/r(ii/6)i 
t  no 

r  ZT  I0^T/2£ui2^e:I?(_^/2i“2J  ~  cos 


(44) 


Ve  note  that  one  half  c f  *-»e  above  formula,  after  also  replacing  { - } 

by  unity,  is  the  generalization  of  the  ray  angle-of-arrival  power  spec  truss 
for  «  «  a  «  L-.  It  yields  the  same  Units  as  (36b)  and  (38b)  with  _  -  0. 

4.2.3  Concluding  Renarks 

The  differences  between  ray  and  interferometer  formulas,  the  latter  given 
by  (38)  and  (40)  and  the  former  by  one  naif  of  the  leading  first  terras  of  each, 
appear  at  first  sight  to  be  unimportant  -  The  curly-bracketed  factors  reduce 
to  unity  when  ' .-/--j-  (°r  in  the  other  case)  exceeds  unity.  Considering 
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fart  e4«  practical  iscerfcrotttcn  to  tot*  $  £.*. ,  5*1  ^ 

*  1  *"**  ffat*  ^«Wac£<.  #  w?  close  to  a.  ^  to*  via*  is  itoefa. 


U  “  **“•  »f  «*«  toarteraror  i.  ****  cm*  b  -  «.  koa*.  «a*  wtor- 
*£^**  [' 11  * mwBlly  terixoattai,  *w  wl(m  dfcas  esc  r bii i  are  <61 

®k: 


£->  7***zeaL  vsSerftarsmetee  {me  %j  nwyiimn,  aai  -  -1>:  15*  aoenre 

rcMBrks  ?cruia  aai  ®§C»>  is  twice  toe  ray  Sff  {*)  for  *11  practical  per- 


tU)  &-'3CTr2cZ  ***ffr*Me*  {me  £*  «*p««t}:  I*  ai,  case  «*  «* 

***  Ci  *  0  for  *11  practical  purposes  to  torato  nrfee  toe  ray  fora*!*  for 

®*(<tt>  *  *  ftiter  }l  -  ^*CC#u/^  ittkcase&t 

fs  not  «U,  aarf  toe  factor  of  (4®>)  for  toe  ©toer,  smell  c*se. 

11  is  tofetoos  toat  *ore  actuate  foraelas  are  retiree  to  toe  u"2/3  re^toe, 
£*=t  toey  caa  ©e  seikei  oet  awertoallj  fr<*  toe  font  *£**»  to  tois  seettoa. 
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